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Climate change and potential selection for 
non-diapausing two-spotted spider mites on strawberry 
in southwestern British Columbia 

D.A. RAWORTH 1 


ABSTRACT 

A validated model of the timing of post-diapause oviposition in the two-spotted spider 
mite, Tetranychus urticae Koch, was used to predict when 50% of strawberry leaflets 
with T. urticae also have T. urticae eggs (IO 0 . 5 ) in each year from 1954 to 2006 at Lang¬ 
ley, British Columbia. This timing was studied in relation to hours of frost occurring 
before and after oviposition. Historically, IO 0.5 occurred before there were frost-free 
days, but there was a clear threshold at 390 h with temperatures < 0 °C after IO 0 . 5 , which 
was not exceeded. This suggests that there is selection pressure for early oviposition, but 
also a limit to the extent of selection. The subzero temperature profile ~1 month before 
oviposition was clearly different from that after IO 0 . 5 . The number of hours with subzero 
temperatures 1 month before oviposition, and the standard deviation of those estimates, 
were negatively correlated with year and indicated that there could be oviposition in 
January - rather than February - by 2015. Cumulative hours with temperatures < 0 °C 
between 27 November (the empirical estimate of the time when T. urticae begins accu¬ 
mulating degree-days for post-diapause oviposition) and 30 April was negatively corre¬ 
lated with year, and extrapolation of a linear regression suggested that there could be 
selection for continuous annual oviposition by 2050. There was considerable variation in 
the data, but considered in combination with published evidence for climate change, 
these results will be important in developing pest management strategies, and further¬ 
more, will impact many aspects of agriculture in the Fraser Valley of British Columbia. 


INTRODUCTION 


Global warming (Intergovernmental 
Panel on Climate Change 2007) will proba¬ 
bly affect arthropod ranges (Gray 2004; 
Logan and Powell 2004; Gutierrez et al. 
2006; Musolin 2007) assuming fixed bio¬ 
logical tolerances to environmental condi¬ 
tions. At the same time, within a home 
range, it will probably also affect the life 
history characteristics of arthropods 
through selection (cf. Bradshaw et al. 
2004). Winter diapause is a key feature of 
temperate arthropods (Danks 2006). Dia¬ 
pause characteristics will probably be af¬ 
fected by global warming, and changes may 
become evident first in areas that have a 
relatively mild but temperate climate, such 
as the Fraser Valley of British Columbia 


(B.C.), Canada. This study considers the 
effects of climate change on the timing of 
initial, post-diapause oviposition (IO) by 
two-spotted spider mites, Tetranychus urti¬ 
cae Koch (Acari: Tetranychidae), on straw¬ 
berry ( Fragaria x ananassa Duch. 
Rosaceae). 

T . urticae females have a facultative 
reproductive diapause. Diapause is induced 
in the pre-imaginal stages, particularly at 
the end of the protonymphal instar, by 
short-day photoperiods (Veerman 1977a). 
Termination is dependent on duration of 
cold rest, temperature, and photoperiod 
during the first few months of diapause 
(Veerman 1977b). During winter, after 
photoperiodic sensitivity is gone, diapause 
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would be sustained and prolonged by low 
temperatures. Under long days at relatively 
high temperatures, diapause can be termi¬ 
nated without a cold rest period (Veerman 
1977b), but under short days a cold rest 
period is required. 

Given a warmer climate there should be 
selection for individuals who reproduce 
continuously, despite short days, and there 
may be selection for individuals who do not 
respond to the initial photoperiodic induc¬ 
tion cues. This idea is supported by the fact 
that there is considerable variation in dia¬ 
pause induction and termination character¬ 
istics among T. urticae strains, probably 
associated with adaptation to local condi¬ 
tions (Takafuji et al. 1991; Koveos et al. 
1999). 

Raworth (2007) developed and validated 
a FORTRAN program (ProglO) that deter¬ 
mined the timing of post-diapause oviposi- 
tion - the day when 50% of strawberry leaf¬ 
lets with T. urticae also had T. urticae eggs 
(IO0.5). The oviposition model for the pro¬ 


gram was based on field samples from 
inland and coastal sites, collected during 
1988 to 2006; estimates of IO 0.5 from 10 
populations during 1988 to 2003 were used 
to calibrate the model, and six independent 
populations during 2004 to 2006 were used 
for validation. 

Here, ProglO and historical weather 
data are used to examine 1 ) where, histori¬ 
cally, T. urticae has placed IO 0.5 with re¬ 
spect to frost, and 2 ) the annual historical 
trends in hours of frost. These trends are 
then used to predict when T. urticae can be 
expected to commence oviposition in Janu¬ 
ary - rather than the current timing in Feb¬ 
ruary - and when T. urticae may go through 
the winter with no reproductive diapause. It 
may be objected that viable leaves must be 
available for feeding and oviposition; how¬ 
ever, in this system mature green leaves 
overwinter, and the mites usually initiate 
feeding and oviposition on these leaves in 
February and March (Raworth 2007). 


MATERIALS AND METHODS 


ProglO and meteorological data were 
used to predict IO 0.5 each year from 1954 to 
2006. ProglO first calculated the tempera¬ 
ture under a strawberry leaf in a commer¬ 
cial field at Langley, B. C. based on hourly 
temperature and cloud opacity data meas¬ 
ured by Environment Canada at Abbots¬ 
ford, B. C. since 1953, and calibration 
equations (Raworth 2007): 

1. Cloud opacity = 1 (full cloud) 

a. 0800-1600 h: y = -8.756 + 1.009 t + 
1.731 /7 - 0.013 Z 2 - 0.076 h 2 + 0.015 t X h 

R 2 = 0.92, P< 0.0001, 425 d.f. 

b. 1700-0700 h: y = -1.104 + 0.864 t - 0.025 
h + 0.010 Z 2 

R 2 = 0.92, P < 0.0001, 862 d.f. 

2. Cloud opacity = 0 (full sun) 

a. 0800-1600 h: y = -54.250 + 1.116 t + 
9.515 h - 0.370 h 2 

R 2 = 0.85, P< 0.0001, 241 d.f. 

b. 1700-0700 h: y = -3.848 + 1.153 t + 
0.066 h + 0.023 i 2 

R 2 = 0.92, P< 0.0001, 554 d.f 
where: y = temperature experienced by the 


mites; t = temperature in a Stevenson 
Screen in the field (t = -1.141 + 1.043 tec , 
where tec = air temperature at the Environ¬ 
ment Canada station); and h = hour (where: 
0800-1600 h = 8 , 9,...16; and 1700-0700 h 
= 17, 18,...24, 25, 26,...31). Temperatures 
under a leaflet in intermediate cloud condi¬ 
tions were determined by linear interpola¬ 
tion between the results provided by the 
equations for full cloud and full sun. The 
timing of IO0.5 was then determined from 
thermal summations > 9.4 °C starting on an 
empirically-derived day, 27 November 
(Raworth 2007), and a thermal requirement 
(y) that was negatively correlated with ac¬ 
cumulated cold-rest hours < 4 °C (x) 
summed from 27 November (Raworth 
2007): 

3. y = 78.3 - 0.0279 x; r 2 = 0.83, P = 0.01, 
4 d.f. 

Equation 3 implies that, for equivalent 
rates of thermal summation, the spider 
mites will stay in diapause longer if they 
have had insufficient cold rest. The timing 
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of IO 0 .5 was compared graphically with the 
cumulative daily hours of frost below 0 °C 
summed between 27 November and 30 
April of the following year. 

To determine what temperature condi¬ 
tions the mites have avoided, annual fre¬ 
quency was plotted against the number of 
hours of frost below 0, -1, -2,...-10 °C that 
remained after IO 0.5 (3-d plot), and the 
hours of frost below 0, -1, -2,...-10 °C that 
occurred one month before mite eggs would 
be observed in the field, 10 January to 10 
February (3-d plot). To determine when 
conditions would be suitable for earlier 
oviposition by T. urticae , the annual sum of 
hours of frost below 0, -2, -4,...-10 °C that 


occurred between 10 January and 10 Febru¬ 
ary were regressed (SAS Institute 2004) 
against year and the regression was ex¬ 
trapolated beyond 2006. Changes in vari¬ 
ability as a function of year were deter¬ 
mined by pooling the latter data for each 
decade and regressing 1 SD for mean hours 
below a given temperature in that decade 
against the median year. Finally, a similar 
technique was used to determine when con¬ 
ditions would be suitable for T. urticae to 
pass through the winter without reproduc¬ 
tive diapause based on the annual sum of 
the hours of frost below 0 °C between 27 
November and 30 April. 


RESULTS AND DISCUSSION 


Cumulative hours of frost below 0 °C 
between 27 November and IO 0 . 5 , from 1954 
to 2006 varied between 400 and 1600 h 
(Fig. IB, triangles). Despite this variation, 
IO 0.5 always occurred before frost-free days 
had begun (Fig. 1A, B), but with a clear 
upper threshold of 390 h of frost remaining 
after IO 0 . 5 . This indicates sensitivity to 
frost, because all the estimates of IO 0.5 are 
clustered below 390 h; it also indicates 
some selection pressure for early oviposi¬ 
tion, because the mites do not wait until 
frost-free conditions occur. As long as a 
female’s progeny can survive and go on to 
reproduce, a female that initiates reproduc¬ 
tion early in the season should have a nu¬ 
merical advantage over one that initiates 
reproduction later. This result should be 
qualified. The 390 h threshold was deter¬ 
mined from a model of post-diapause ovi¬ 
position in southwestern B. C.; it would not 
be expected to apply to T. urticae popula¬ 
tions that are adapted to different local con¬ 
ditions in other temperate regions. Such a 
generalization would require further re¬ 
search. 

Comparison of the subzero temperature 
profiles before and after IO0.5 suggests that 
T. urticae in southwestern B. C. is able to 
tolerate >200 h at -2 °C, but only 20 h at -6 
°C (Fig. 2). Lower temperatures for longer 
periods (Fig. 3) have been avoided. How¬ 


ever, the number of hours with tempera¬ 
tures below 0, -2, ... -10 °C, 1 month before 
T. urticae normally commences oviposi¬ 
tion, have declined significantly since 1954 
(Fig. 4). Taking -8 °C as a critical tempera¬ 
ture associated with no oviposition (Fig. 2), 
the data indicate that an average year will 
have zero hours <-8 °C by 2015 (Fig. 4). At 
this time there could be reduced selection 
pressure against emergence, and hence ovi¬ 
position in January rather than February. 
Because the regression predicts hours at a 
given temperature in an average year, one 
would expect some hours at temperatures 
<-8 °C in some years, and possible mortal¬ 
ity of early-emerging mites. However, the 
variation in subzero temperatures during 
this 1 month period has also declined since 
1954 (Fig. 5) so that there will be less un¬ 
certainty about subzero conditions in 2015 
than there was in 1957, and reduced selec¬ 
tion pressure against early emergence. The 
regression of total hours of frost below 0 °C 
from 27 November to 30 April against year 
suggests that the number of hours of frost 
will decrease to the threshold of 390 h by 
-2050 (Fig. 6 ). At this point, selection for 
early oviposition could result in reproduc¬ 
tion by some individuals right through the 
winter in southwestern B.C. 

The objective of this study was to exam¬ 
ine general patterns in the initiation of post- 
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Figure 1 . Cumulative hours with temperatures < 0 °C under a strawberry leaflet in a commer¬ 
cial field from: (A) start date to 30 April versus start date (dots); (B) 27 November to the date 
when 50% of strawberry leaflets with T. urticae also have T. urticae eggs (IO 0 . 5 ) versus IO 0.5 
(triangles); and (B) IO 0.5 to 30 April versus IO 0.5 (circles), each year from 1954 to 2006. The 
horizontal line represents a threshold number of hours < 0 °C remaining after IO 0 . 5 . 
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Figure 2 . Number of years with a given number of hours at subzero temperatures under a 
strawberry leaflet in a commercial field, after the date when 50% of strawberry leaflets with T. 
urticae also have T. urticae eggs. 
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Figure 3. Number of years with a given number of hours at subzero temperatures under a 
strawberry leaflet in a commercial field, ~1 month before the date when T. urticae commences 
oviposition (10 January to 10 February). 



Figure 4. Hours with subzero temperatures under a strawberry leaflet in a commercial field, 
between 10 January and 10 February, versus year. Solid lines are linear regressions (r 2 = 0.10, 
0.08, 0.08, 0.09, 0.10, 0.14 and P = 0.02, 0.03, 0.03, 0.02, 0.02, 0.004) for temperatures < -10, 
-8, ... 0, respectively; data for temperatures <-4 °C are shown. 
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Figure 5. One SD of mean hours with subzero temperatures under a strawberry leaflet in a 
commercial field, between 10 January and 10 February, versus median year. Solid lines are 
linear regressions (r 2 = 0.95, 0.92, 0.89, 0.87, 0.93, 0.80 and P = 0.0008, 0.002, 0.005, 0.007, 
0.002, 0.02) for temperatures < -10, -8, ... 0, respectively. 



Figure 6. Hours with temperatures < 0 °C under a strawberry leaflet in a commercial field, 
from 27 November to 30 April, versus year. Solid line is a linear regression y = 19600. - 9.37 
x; r 2 = 0.22; P < 0.001. 


diapause oviposition in T. urticae in rela¬ 
tion to field temperatures. A number of 
thresholds were observed for spider mite 


populations in southwestern B. C.: a maxi¬ 
mum of 390 h of frost after IO 0 . 5 ; post-IO 0 .5 
tolerance to >200 h at -2 °C, but only 20 h 
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at -6 °C; and -8 °C was taken as the critical 
temperature after IO 0 .5 associated with no 
oviposition. These thresholds could be con¬ 
firmed by further work, but this was not the 
point of the study; regardless of the exact 
value of the thresholds, the patterns are 
clear. Oviposition is initiated before frost- 
free days occur - in ProglO and ipso facto 
in the field data used to build ProglO; the 
number of hours of frost 1 month before 
oviposition, and variation in those estimates 
among years, has decreased significantly 
during the last half century; and the number 
of hours of frost from 27 November 
through 30 April have also decreased sig¬ 
nificantly over the years. 

It is clearly risky to extrapolate from a 
linear regression based on data with consid¬ 
erable scatter (Figs. 4, 6); the relationship 
may be negative but asymptotic at, for ex¬ 
ample, 600 h frost (Fig. 6). However, there 
is additional evidence. Although it is diffi¬ 
cult to attribute observed temperature 
changes to natural or human causes at 
smaller than continental scales because 
factors such as land use change and pollu¬ 
tion complicate the picture (Inter¬ 
governmental Panel on Climate Change 
2007), the trends observed in the current 
study are consistent with the global warm¬ 
ing scenario. The mechanism driving global 
warming, namely increasing levels of 
greenhouse gases (GHG) (N 2 0, CH 4 , and 
C0 2 ) is well established, and ‘With current 
climate change mitigation policies and re¬ 


lated sustainable development practices, 
global GHG emissions will continue to 
grow over the next few decades’ (Inter¬ 
governmental Panel on Climate Change 
2007); a global temperature change of +0.2 
°C per decade is projected. Therefore, the 
linear extrapolations in Figs. 4 and 6 may 
be reasonable. Time will tell, however if 
correct, oviposition in January should be 
observable within the next decade. 

Despite the many uncertainties in this 
study, there is sufficient evidence for earlier 
post-diapause oviposition in T. urticae 
within the relatively near future in south¬ 
western B. C. to at least consider it in plan¬ 
ning spider mite monitoring and manage¬ 
ment activities. This would evolve naturally 
into planning for continuous annual ovi¬ 
position should that occur in 4 to 5 decades. 
Furthermore, a relatively rapid reduction in 
the number of hours of subzero tempera¬ 
tures during the winter will have significant 
implications for many aspects of agriculture 
in the Fraser Valley, including arthropod 
pest and disease management, crop produc¬ 
tion, and crop selection in both field and 
greenhouse environments. These effects 
need to be considered carefully by growers, 
pest managers, researchers, and government 
planners at Provincial and Federal levels 
with studies and approaches to address po¬ 
tential problems evolving as trends become 
increasingly certain. 
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